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The Rayleigh distribution is the distribution of the sum of a large number of coplanar 
(or time) vectors with random amplitudes and uniformly distributed phases. As such, it is 
the limiting case of distributions associated with more general vector sums that arise in 
practical problems. Such cases arc the following: (a) The phase distributions of the vector 
terms are not uniform, e.g., in the case, of scattering from rough surfaces; (b) One or more 
vector terms predominate, their mean square value not being uegligible compared to the 
mean square value of the sum, e.g., in the case of signals propagated in cities, meteor-scatter, 
and atmospheric noise; (c) The number of vector terms is small, e.g., in radar returns from 
several close targets; (d) Pie number of vector terms is itself random, e.g., in atmospheric 
turbulence, meteor-scatter and atmospheric noise. The resulting distributions for these 
cases and their deviations from the Rayleigh distribution will be considered. 



1. Introduction 

In many problems arising in radio wave propaga- 
tion the resultant field is formed by the super- 
position or interference of a number of elementary 
waves: 



Ee id =J^E J e i ^, 



(1) 



where the K } and the </> ; and even n may be random, 
and the distributions of the individual terms of the 
sum need not be identical. One is then faced with 
the problem of determining the distribution of E 
(and sometimes of 0) if the distributions of the E j7 
the 4>j and n are known. It will be assumed that 
the terms of the sum (1) are mutually independent. 
The sum (1) may also be regarded as the sum of 
coplanar vectors. 

In its most elementary form, when n is a large 
constant, the E i are all equal to the same constant, 
and the </> ; are all uniformly distributed from to 
2w, the problem was solved by Rayleigh [1896] and 
leads to the well-known Rayleigh distribution 






(2) 



where x = (E 2 )=nEj and the brackets ( ) denote the 
mean value; or 
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P(E/E JlMS >R) = e- R2 ; E RMa = Vi- 



(3) 



The present paper considers (1) under more 
general conditions. In principle the problem can 
always be solved by resolving (1) into its rectangular 
components (i.e., real and imaginary parts) x and y y 



x=E cos 6=^2 Hj cos <t>j='^2^i 

J=U " .7=1 

y=E sin 6= ... 



(4) 



finding the joint probability density W(x,y), and 
re transforming to polar coordinates through 

p(E) =E \W(E cos 6, E sin d)dQ 

Jo 

/» 00 

EW(E cos 6,E sin 6)dE. (5) 



The various cases that may arise are conveniently 
classified according to whether or not the sums (4) 
satisfy the Central Limit Theorem, i.e., whether or 
not x and y are normally distributed as n— >°°. 
Each of these two cases again includes several 
further possibilities. 

Let D(x) denote the variance of x; then the distri- 
bution of x will tend to a normal distribution as 
n— >oo provided that the Xj are mutually independent 
and that 



D( Xj ) 



lim : n 



= for all j, 



(6) 
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with a similar statement for y. (This statement 
usually suffices for engineering purposes; for a more 
rigorous enunciation of the Central Limit Theorem 
and Lindeberg conditions cf. Gnedenko and Kol- 
mogorov [1954].) Condition (6) essentially means 
that none of the terms x } must predominate in the 
resulting sum x. However, if the Xj are themselves 
normally distributed, then x will, of course, be also 
normally distributed even if (6) does not hold. 

In most (but not all) applications the <j>j are 
distributed uniformly between and 2t or in an 
equivalent manner; i.e., the phase distribution 
w<t,((j)) is such that 



S ^(*+2fcr)=^+e(0); |e(4>)|« 



2tt 



(7) 



A vector with such a Uniformly Distributed Phase 
will be called a UDP vector. The sum of UDP vec- 
tors is obviously itself a UDP vector. If the terms 
in (1) are UDP vectors, then 



(E 2 )=(12 Eje*** S E k e-*+* 



=I2(E>)+T J {E j E k )(e^-^). (8) 

But for j^k, 

( e <(*,-**>)=_ €* ( ^-** ) ^d**=0. (9) 

47T" Jo JO 

Substituting (9) in (8) we find the important 
relation 






(10) 



valid for UDP vectors regardless of the value of n 
or the distributions of the Ej (possibly all different). 

2. Rayleigh Distribution 

If the terms in (1) are UDP vectors, then from (4) 
<*>=<y>=0 (11) 

Z?(a!,)=I?(y,)=^J o ) o EUos'hdEjd^iEf). 

(12) 
Then using (10), condition (6) becomes 



lim } E ^ =0 for all ;'. 



(13) 



If (13) is satisfied and n is large, x and y can be 
approximated by a normal distribution with mean 
value zero and the same variance. The integration 
(5) then leads to the Rayleigh distribution (2) with 



7 = 1 



(14) 



Thus a Rayleigh vector is a UDP vector whose x 
and y components are distributed normally with 
(x) = (y) = and D(x) =D(y)=s/2. From _ this it fol- 
lows that the sum of any number of Rayleigh vectors 
is itself a Rayleigh vector. 

A Rayleigh distribution will thus be found when- 
ever the resultant field is composed of a large number 
of UDP vectors and (13) is satisfied. 



3. Nonuniform Phase Distributions 

In a number of applications the phases <f>j in (1) 
are not distributed uniformly as in (7) , but fluctuate 
about some privileged value. This will occur in 
scattering from rough surfaces (e.g., rough layers 
in the atmosphere) for small roughness or small 
grazing angles. Since the terms in (1) are now not 
UDP vectors, (11) and (12) will not hold. How- 
ever, if (6) holds; x and y will still be normally dis- 
tributed. If the phase distributions are symmetrical 
about zero, then (y) = 0. By the usual rules of 
probability theory one then finds the quantities 

a=<3>> 8i=D(x), s 2 =D(y) 
and the integration (5) yields 

Vsi«2 L 2s, 4.s 1 s 2 J 



TO = 



■'-(S'' s Mr) < 15 > 



where I m is the modified Bessel function of order 
m and e =l, e m =2 for m^O. Details of the pro- 
cedure and curves of (15) will be found in [Beckmann, 
1962a]. 

The general distribution (15) simplifies in certain 
special cases. If a = 0, but s^s, then (15) reduces to 



p(E) = -=exip 



r si+s 2 t^i T A2--S1 

L 5* J His; 



f z 



(16) 



a distribution derived directly by Hoyt [1947]. 

On the other hand, if s 1 =s 2 =tzS, but a^O, then 
(15) reduces to the Nakagami-Rice distribution 

a 2 +E> 



fin 2E 
p(E) =— exp 



i-m< 



2aE\ 



(17) 



a distribution derived by Rice [1944 and 1945] 2 
and further analyzed, e.g., by Norton et al. [1955], 
and Zuhrt [1957]. The distribution (17) is obtained 
when a constant vector (Ei=a } 0i = O) is added to a 



2 The distribution was originally derived by Nakagami in 1940. A summary 
and bibliography of the work on this and related topics by Nakagami and other 
Japanese scientists will be found in [Nakagami, I960]. 
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Rayleigh vector, for the x and y components of the 
sum will then obviously be distributed normally 
with (x)=a, {y) = 0, D(x)=D(y)=s/2 (where s is 
the mean square value of the Rayleigh vector) 
just as assumed in deriving (17). It should be noted 
that (17) will equally well hold for E } exp (ifa) a 
UDP vector with constant amplitude E l =a; this 
may be shown by measuring the p liases from fa as 
a reference phase: the distributions of ^ = 4,.— fa 
will for j 5^1 remain uniform as before, whereas </>i = 0, 
thus reducing to the same conditions under which 
(17) was derived. 

The Rayleigh distribution, as may easily be 
verified, is the limiting distribution of (15), and of 
its special cases (16) and (17), fora=0; s l =s 2 =s/2. 

If the Central Limit Theorem may be applied to 
(4), so that x and y will be distributed normally, 
then in the most general case (corresponding to 
asymmetrical phase distributions) we have four 
parameters: 

(z)=a; (y)=Pi D(x)=su D(y)=s 2 . (18) 

The integration (5) then leads to 

ivi p — £ CO 



COS 



[~2m (arc tan ^Yl (19) 



where 



Si 



4:SiS 2 2si 2s 2 ' 4.v r vo 



,2. 



To. 



iQ=^ B 



Si 



$2 

(20) 



The distribution (19) was found by Nakagami 
[I960]; it is the most general distribution for the case 
when the Central Limit Theorem is applicable to 
(4); for |S=0 and hence R=0, it reduces to (15). 
The Rayleigh distribution is again obtained from 
(19) for a=P=0; s } = s 2 =s/2. 

If the Central Limit Theorem is not applicable to 
(4), this may be for one of the following reasons: (a) 
condition (6) is not satisfied (this will be considered 
in sees. 4 and 5); (b) the number of terms n in (1) is 
not large (sec. 6) or random (sec. 7). 

4. Dominant Terms 

For UDP vectors, which we shall henceforth as- 
sume, (6) reduces to (13). If the number of inter- 
fering waves n is large, but finite, we may replace 
(13) by 



<^S>«S <JSJ> for any j. 



(21) 



Now if one (or more) of the interfering waves is 
powerful, so that its power is not negligible when 
compared to the total power, (21) will be violated 



for one (or more) j. If more than one of such waves 
are present, we may sum them by standard methods 
(convolutions, characteristic functions) and regard 
this partial sum as one wave. Since assumption (21) 
will hold for all other j, the remaining terms will add 
up to a Rayleigh vector, so that the problem reduces 
to finding the distribution of the sum of a UDP 
vector with random amplitude E x and a Rayleigh 
vector. 

This problem may be solved directly from first 
principles by (4) to (5) or more quickly by random- 
izing a=Ei in (17) and using the theorem of total 
probability: if the density of a is w(a), then (17) 
gives the density p(E\a), so that the required total 
probability density is 



p(E) = 



2E 



/; 



w(a) exp ( - 
. 1 \ 



cr+E 2 \ r /2aE\ 



> 



-jda. 



(22) 



It may be verified from (22) that p{E) will 
approach a Rayleigh distribution for(a 2 )<C < 0, as was 
to be expected. The complement of the distri- 
bution function of (22) is 

P(E>R) = f " p(E)dE= r w(a)j(R, a)da (23) 

J R Jo 

where the order of integration has been reversed with 

(24) 

Now if It is large (R>y>s/a), the Bessel function 
in (24) may be replaced by its asymptotic expression: 
a saddle-point integration then leads to 



AB.4-i['-(^)} 



(25) 



Now (25) changes its value from <().()! to >0.99 
near the point a =R within an interval Aa=2.3V^R, 
tending to zero below and to unity above that 
interval; for fi» Vs we may therefore well approx- 
imate (25) by 



/(£,«): 



for a<R 

1 for a >R. 



(26) 



Substituting this value in (23) we find 



P(E>R) « (" w(a)da=P(a>R) (for B»Vs). 

J R 

(27) 

Hence for /?>>->/<$ the distribution of a random 
vector plus a Rayleigh vector will approach that of 
the random vector alone. This effect may be 
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observed in several cases in radio wave propagation. 
One of these is the field strength of VHF and UHF 
in cities and other built-up areas, where the total 
signal may be due to a direct wave (attenuated at 
random as it is transmitted through walls and other 
building materials) onto which large numbers of 
reflected waves are superimposed. The resulting 
amplitude (which is constant in time, but random 
when measured at different places) is then dis- 
tributed as in (22). A survey conducted at various 
parts of the city of Prague showed that in most 
areas w(a), the distribution of the attenuated 
direct wave, is lognormal; the same result need not 
necessarily hold for cities differing in character from 
the above (brick or concrete houses, five to six 
stories high, streets relatively narrow and not 
forming a regular pattern). By analyzing the 
resulting distribution the propagation mechanism 
may thus be investigated (separation of reflections 
and attenuation). 

5. Converging Variances 

In at least two cases met in propagation theory, 
meteoric forward-scatter and noise due to atmos- 
pherics, the signals arrive at random intervals of 
time with a random amplitude which then decays 
exponentially. The signals are mutually independ- 
ent and their phase makes them UDP vectors. 
Since an exponentially decaying signal never vanishes 
completely, there is an infinity of residual signals 
present at any time; but since the power at any time 
is finite, the infinite series of signals must converge 
and the denominator of (13) will not tend to infinity 
with n (this can also be shown mathematically). 
Hence (13) will not be satisfied for any j, the Central 
Limit Theorem cannot be used, and the resulting 
amplitude distribution cannot be a pure Rayleigh 
distribution. 

To solve the problem rigorously one therefore 
has to return to first principles, e.g., by finding x 
and y in (4) through their characteristic functions 

X x iv) =Xy (v) =IL t^V, (14>\ dEjw(Ej)e ivE J C08 « 
j=i (2t) j Jo Jo 



co n<* 

= n 

j=i Jo 



wiEJJoivEddE, 



(28) 



where Ej is the amplitude of the^'th decaying signal, 
which depends on two random quantities: the time 
tj elapsed since the signal attained its peak value, 
and that peak value E p : 



E^E v e- 1 ^ 



(29) 



with a the time constant of decay. Since the number 
of signals per unit time is Poisson-distributed (about 
an average N), t x has an exponential distribution 



w 1 (t 1 )=Ne~ 



(30) 



and the distribution of tj is given by a^'-fold convolu- 
tion of (30), which leads to 






(31) 



The probability density of 7?=exp (—tjja) is then 
found from (31) by a simple transformation: 

Wj J {j-l)\B aN + l ( 1 <5<oo). (32) 

If the density of E p is \(E P ), then the distribution 
of Ej is found from (30) as that of a product of two 
random variables: 

w(Ej) =f~EMB j E j )W j (B j )dB J . (33) 

Substituting (33) in (28) will in general lead to 
great computational difficulties, which may be 
overcome by the following approximation. 

From (33) we find 

(Ef)=^E 2 j w{E j )dE 1 (34) 

(35) 



and 






This series will converge and hence violate (13); 
however, if the series converges quickly (Na<C<Cl), 
then (21) will still hold for jj*l; if this term is 
excluded, the rest of the series may thus be approxi- 
mated by a Rayleigh vector. Physically this means 
that if the time constant of decay is sufficiently 
short for the signal to decay to a low value before 
the next signal arrives (average interval is 1/N), 
then the total signal at any time will be dominated 
by the last signal (or possibly the last few signals), 
whilst the remnants of all previous signals will 
combine to form a low-power Rayleigh vector. 

Tims under these circumstances this case may 
be reduced approximately to the one in section 4; 
the required distribution is thus given by (22) with 
w(a)=w(E 1 ) given by (33) and s=(E 2 )-(Ef) found 
from (34) and (35). 

The distribution \(E P ) is given by the physical 
nature of the problem. For atmospherics, \{E P ) 
may be shown to be lognormal [Beckmann, 1964]; 
in spite of the several approximations involved, 
the agreement of the distribution as derived above 
and the measured distribution is very good as shown 
by figure 1. An analysis of the distribution permits 
the effect of propagation conditions to be separated 
from that of lightning activity in the total random 
atmospheric noise. 

6. Small Number of Independent 
Components 

If the number n of independent interfering waves 
is small (e.g., the total radar signal returned from a 
small number of independent targets in the same 
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Figure 1. — Amplitude-probability distribution of atmospheric 
radio noise. 

Circles: values measured by Crichlow et al. [19601, on 13.3 kc/s at Boulder, 
Colo., October 6, 1958. Full Curve: distribution computed as in (33) through 
(35) for <7=.1.67, iVc=0.01. Broken curve interpolated. (Cf. Beckmann [1904].) 

area), the Central Limit Theorem cannot be applied 
and little can be said about the distributions of ./■ 
and y in (4) beyond the statement that they have 
to be determined from the Wj(Ej) by convolutions 
or characteristic functions. However, we may ask 
the opposite question: how large must n be in (1) 
for the Rayleigh distribution to be a good approxima- 
tion for p(E) in engineering practice? The answer 
will obviously depend on the distributions Wj(Ej), 
which we here assume all identical and equal to 
w(Ej). In that case one may derive the formula 
[Watson, 1944; Levin, 1960] 

V {E)=eT ( I ^ wiE^J.iuE^dEXuJoiE^du. 

(36) 

Expanding J^E/u) in a series and integrating 
term by term, rearranging in ascending powers of 
(l/n) and using the 2d and 4th initial moments of 
w(Ej), i.e., 



m 2 = <JBJ> = 



r 



EfoiEddE,; 



f 



Eju)(Ej)dEj (37) 



we find an expression which on integration over E 
yields [Levin, 1960, pp. 184-187], 

(38) 

For w— >oo this leaves a pure Rayleigh distribution; 
for finite n, the Rayleigh distribution will be a good 
approximation if the second term in the square 
bracket will be small compared with unity; i.e., the 
required criterion is 



n» T 



0-2 1 



(39) 



where p^mJm^iEfiKE*) 2 . 

7. Random Number of Terms 

In most cases met in wave propagation through 
random media the number of interfering waves n is 
not constant; the number of scatterors such as 
turbulence cells in the atmosphere changes from 
moment to moment; the number of effective reflec- 
tors in cities or irregular terrain changes from loca- 
tion to location, etc. Thus n itself will be random; 
if its distribution (discrete for integers only) is 
/'(//), then the distribution of x and y in (4) will be 

p x (x) = S P(n) S Ej cos (/>,, (40) 

3t=0 j=l 

Now if P(n) assumes appreciable values only for 
large n, the terms of the w-sum will be normal, hence 
x will be normal and E will be Rayleigh-distributed. 
A more detailed investigation [Beckmann, 1962b] 
shows that the distribution of Ewill closely approach 
a Rayleigh distribution as the condition 



<«> 2 " 



>i 



(41) 



is more nearly satisfied; it is also shown that for 
a given distribution P(n) the deviation of P{E) 
from a Rayleigh distribution will always be greater 
for large E than for small E. 

As a rule n is Poisson-distributed about its mean 
value (n), i.e., 



P(f0=^6<-»>. 



From (42) we find (n 2 ) 
(n 2 ) 



{nf (n) +1 



'-{n) + {n) 2 ; hence 



1 



(42) 



(43) 



which will approach unity as required by (41) for 
(n)»l. Thus if P(n) is given by (42), a large 
mean value is sufficient to make the Rayleigh dis- 
tribution a good approximation for p(E). It should 
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be noted that for large (n), (42) is equivalent to a 
normal distribution with D(n)=(n) [Levin, I960]. 

For a general distribution P(n), however, a large 
mean value (n) is not sufficient to guarantee 
that p(E) in (1) and (4) will tend to the Rayleigh 
distribution, but the criterion (41) must be satisfied. 

8. Distribution of the Resulting Phase 

The distribution of 6 in (1) is found from the second 
relation in (5). For a sum of UDP vectors, the 
resulting phase will of course again be uniformly 
distributed over an interval of length 2tt; for non- 
uniform phase distributions of the vector terms we 
introduce the quantities 



B- 



and then obtain from (5) 



;K= 



V's 2 
ft 



(44) 



P»(fi> 



j^ e -iB 2 a+K 2 ) 



2t(K 2 cos 2 0+sin 2 0) 



(45) 



where 



G=BK cos B 



4 



l+K 2 



2{K 2 cos 2 0+sm 2 6) 



(46) 
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